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ABSTRACT 



o 



. Context. The 3rd year WMAP results mark a shift in best fit values of cosmological parameters compared to the 1st year data and the concor- 
• dance cosmological model. 

Aims. We test the consistency of the new results with previous constraints on cosmological parameters from the HIFLUGCS galaxy cluster 
f*^) sample and the impact of this shift on the X-ray luminosity-gravitational mass relation. 

' Methods. The measured X-ray luminosity function combined with the observed luminosity-mass relation are compared to mass functions 
CZ^ predicted for given cosmological parameter values. 

Results. The luminosity function and luminosity-mass relation derived previously from HIFLUGCS are in perfect agreement with mass 
functions predicted using the best fit parameter values from the 3rd year WMAP data (Q. m = 0.238, <x 8 = 0.74) and inconsistent with the 
, concordance cosmological model (fl m = 0.3, erg = 0.9), assuming a flat Universe. Trying to force consistency with the concordance model 
O ^ requires artificially decreasing the normalization of the luminosity-mass relation by a factor of 2. 

Conclusions. The shift in best fit values for Q m and cr 8 has a significant impact on predictions of cluster abundances. The new WMAP results 
are now in perfect agreement with previous results on the Q m -cr 8 relation determined from the mass function of HIFLUGCS clusters and other 
X-ray cluster samples (the "low cluster normalization"). We conclude that - unless the true values of fi m and cr 8 differ significantly from the 
£-3 3rd year WMAP results - the luminosity-mass relation is well described by their previous determination from X-ray observations of clusters, 
• • ■ with a conservative upper limit on the bias factor of 1.5. These conclusions are currently being tested directly in a complete follow-up program 
^ ' of all HIFLUGCS clusters with Chandra and XMM-Newton. 
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1. Introduction iReiprich et al.ll2004T) . gravitation al lensing masses may be af- 
fected by projection effects (e.g.. lClowe et al.l2 004). In the era 

The abundance of local galaxy clusters is very sensitive to the of precision cosmo i og y an d with the realization of very tight 

matter density, Q m , and the amplitude of density fluctuations, constraints on cosmo logical parameters especia lly from cos- 

expressed as cr 8 , in the Universe. This has long been known and mic microwave background measurements (e.g.. lSpergel et"aH 

the method has been applied often to constrain cosmological | 2 Q03|) , these difficulties have lead to the idea to turn the argu- 

parameters (e.g,,|Henry & Arnaud 1 99 1 ; Bahcall & Cen 1 992j ment and assume the cosmo i ogica i mode l to be known 

[19931 | Pierpaoli et al. || 2003l | Voevodkin & Vikhlinin| | 2004 | ) . In precise iy an d use the abundance of clusters to study cluster 

practice, this test requires detailed knowledge about the num- physics . For example, the X-ray luminosity-gravitational mass 

ber density of clusters and their gravitational mass. The number (Lx _ Mtot) relation may be derived by combining the measured 

density as a function of, e.g., X-ray luminosity can be measured x . ray luminosity function with the assumed cluster mass func- 

relatively easily while the accurate determination of the cluster t j on 
mass is difficult and subject to systematic errors. For instance, 

mass determinations based on galaxy velocity dispersion and In the last years, a broad consensus has developed that we 

X-ray emission of the intracluster gas m ay be affected if the live in a concordance Universe (Q m = 0.3, Q A = 0.7, <x 8 = 

cluster is undergoing a merger event (e.g jRandall et alJl2002l 0.9), based on consistent measurements from a multitude of 

completely independent observations (from cosmic microwave 
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* Founded by merging of the Institut fur Astrophysik galaxy clusters etc., e.gjBahcall et al.ll999UWang et al.l2002t 
und Extraterrestrische Forschung, the Sternwarte, and the lTurned l20 02albl) : althou gh the cr s constra int remained some- 
what controversial (e.g., Wan g et al .120031) . Then, the case was 
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Fig. 1. Galaxy cluster mass functions. Filled circles: 
HIFLUCCS data (63 clusters); dashed line: best fit to 
HIFLUGCS data; upper solid line: predicted by concor- 
dance model (Q m = 0.3, <x 8 = 0.9), which was a very good 
fit to the WMAP 1st year data; lower solid line: predicted by 
best fit to the WMAP 3rd year data (Q m = 0.234, cr 8 = 0.74). 
The concordance model predicts too many (>200) clusters in 
the survey volume while the n ew cosmo l ogy pr edicts about the 
correct number. Adapted from lReiprichl j200a. Fig. 10.87). 

further strengthend when the 1st year results from WMAP re- 
turned more evidence for a flat Universe for reasonable values 
of the Hubble constant, and best fit values that are very close 
to the concordance values (Q m = 0.29, cr 8 = 0.9. ISpergel et alJ 
120031 their Tab. 2). 

Although the best fit Q m and <x 8 values changed only within 
the uncertainties between the WMAP 1st year and 3rd year 
data, this has nonetheless dramatic consequences for predic- 
tions of cluster abundances, since these are very sensitive to 
even small changes in these values. This is illustrated in Fig.^ 
where the cluster mass functions for the concordance cosmol- 
ogy and the new b est fit values (O m = 0.238, <x 8 = 0.74, 
Sne rgel et alJ l2006t from now on referred to as "new cos- 
mology") are compared. Also shown are data points from the 
HIFLUGCS mass function, determined with individual mass 
measurements for 63 galaxy clusters, assuming the intracluster 
gas to be in hydrostatic equilibrium with the gravitational po- 
tential ( Reiprich & B6hringerl l2002). The concordance model 
predicts more than 200 galaxy clusters in the survey volume 
while the new cosmology predicts ~50 clusters. It is to be ex- 
pected that this change in best fit cosmology has significant 
consequences when deriving the cluster Lx-M tot relation by 
combining the measured luminosity function with the predicted 
mass function. The aim of this Letter is to illustrate these con- 
sequences and suggest a current best estimate of the Lx-M tot 
relation. 

2. Measuring and Predicting the Cluster L x -M tot 
Relation 

Th e observed cluster Lx-M tm rel ation was first studied in detail 
by iReiprich & BohringeJ ||2P02,RB02) using the X-ray flux- 
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Fig. 2. The best fit model to the WMAP 3rd year data 
lies exactly at the intersection of the Q m -cr 8 relations from 
HIFLU GCS and COBE. COBE relation from iBunn & White! 
lll997l) . calculated using a code provided by T. Kitayama. 
Adapted from Reiprich (2003). 



limited HIFLUGCS sample. It is difficult to accurat ely simu- 
late th e X-ray luminosity of galaxy clusters (e.g. iFrenk et all 
1999) but it is reassuring that high resolution cosmological 
hydrodynamical simulations that include cooling and non- 
gravitational heati ng are able to r eproduce this observed Lx- 
M tot relation (e.g. jTornatore et alJ2003h . 

Since the scatter in this re l ation is finite, RB02 
(Section 4 2), as well as [Reiprich d200ll Section 6.4) and 
llkebe et all |?002, Section 2.3) mention explicitly that the nor- 
malization might be slightly biased high because, in a flux- 
limited sample, clusters with higher luminosity for a given 
mass have a higher likeliho od to be inc l uded in the sample. 
In a recent detailed study, IStanek et alJ J2006I) quantify this 
bias by predicting the Lx-M tot relation using cluster number 
counts and the scatter in the luminosity-temperature relation, 
assuming the cluster mass function in a concordance Universe 
(Cl m = 0.3, Qa = 0.7, cr 8 = 0.9). They find the normalization 
measured by RB02 is biased high by a factor of 2 - a significant 
bias. 

The mass function measured with HIFLUGCS is not con- 
sistent with a concordance Universe (Fig.[Q. RB02 found the 
relation <x 8 = 0.43O m 38 , which surprisingly implied <r 8 = 0.68 
for a fiducial O m = 0.3, and discussed this "low cluster nor- 
malization" for the first time in detail. This relation is shown 
in Fig. El Also shown are the results from the COBE 4 year 
data iBunn & White] 1 19971) and the WMAP 3rd year best fit 
value. The new WMAP result lies exactly at the intersection of 
the COBE and HIFLUGCS results. The fact that the WMAP 
value lies exactly on the HIFLUGCS curve also implies that 
the new cosmology predicts the correct cluster abundance. This 
suggests that the normalization of the Lx-M tot relation mea- 
sured for the HIFLUGCS clusters may actually not be biased 
very strongly. 

To test this quantitatively, we determined "quasi-mass func- 
tions" from the HIFLUGCS luminosity function. In de- 
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Fig. 3. Quasi-mass functions estimated from luminosity func- 
tion and luminosity-mass relations. Filled circles: Lx-M tot re- 
lation as measured with HIFLUGCS; open squares: normal- 
ization of this relation artificially divided by a factor of 2. 
The solid lines indicate the respective best fit mass functions: 
(Q m = 0.22, erg = 0.74) and (Q m = 0.29, cr 8 = 0.88); the 
former result is very close to the new cosmology favored by 
WMAP and the latter very close to the concordance model. 



tail, the availability of individually measured masses for all 
HIFLUGCS clusters was ignored and the mass for each clus- 
ter was instead estimated from its luminosity and the Lx-M tot 
relation (M\L in Tab. 8 of RB02). With these "masses," the 
quasi-mass function was constructed. Then, predicted mass 
functions were fit to the quasi-mass function, following the pro- 
cedure described in detail in RB02, resulting in best fit values 1 
Q m = 0.22 and <x 8 = 0.74 - almost p recisely the values fa vored 
by the new cosmology (Fig. [5}- Since Stane k et alJj200 6) argue 
that the normalization of the HIFLUGCS Lx-M lot relation is 
overestimated by a factor of 2, we divided the normalization 
by 2 and repeated the above procedure. The new best fit values 
were Q m = 0.29 and <x 8 = 0.8 8; i.e., almost exactly the concor- 
dance cosmology assumed bv lStanek et alJ J200cj) . Therefore, 
assuming the new cosmology from the WMAP 3 rd year d ata, 
we c onclude that the factor of 2 bias found by IStanek et alJ 
(2006) is due to assuming the "wrong" cosmological model. 
This conclusion is qualitatively consistent with their relation 
Lisfl °c cr g 4 (note that the "In" in the relation they give in their 
abstract and in their derivation in Section 4.2 seems to be a 
typo), which implies a factor of 3 reduction in normalization 
when going from cr 8 = 0.68 to <x 8 = 0.9. Note that in a revised 
version of their paper, Stan ek etaD will offer a "compromise" 
solution, assuming <x 8 = 0.85 (for Q m = 0.24), resulting in 
moderate scatter in the Lx-M tot relation and a smaller bias (A. 
Evrard, priv. comm.). 

1 Note that llkebe et aljl2002l) and RB02 demonstrated that even for 
a local cluster sample; i.e., without information on cluster abundance 
evolution, the degeneracy between Q. m and crg can be broken through 
the shape of the mass function (the resulting direct best fit values from 
the RB02 mass function are consistent within the statistical uncertain- 
ties with the new cosmology). 



The absence of a strong bias found here also implies that 
the intrinsic scatter in the Lx-M tot relation may indeed be rel- 
atively small. Assume, unrealistically, that there is zero intrin- 
sic scatter. This would imply a bias factor of 1 (i.e., no bias 
at all) for cluster samples drawn from a given mass distribu- 
tion, regardless of how the sample was selected, because all 
clusters with a given mass have exactly the same luminosity. 
The larger the scatter the larger the bias in flux- limited sam- 
ples, t herefore, a small bias indicates small scatter. lzhang et all 
(2006) use 14 cluste rs from an almost volume-complete sam- 
ple. As mentioned in Reiprich (2001), in truly volume-limited 
samples there should be no bias in the normalization of the 
Lx-Mtot relation, regardless o f the amount of intrinsic scatter. 
In their XMM-Newton analysis Z hang et al.l confirm normaliza- 
tion and scatter of the RB02 luminosity-mass relation, corrob- 
orating our result here that the bias is small. 

Note that we do not argue here that the bias factor in the 
normalization of the Lx-M tot relation estimated from flux- 
limited samples is exactly 1 (see our re ferences above), but 
rather that a bias factor of 2 as found bv lStaneket ail (2006) 
is a strong overestimate unless the new WMAP results are 
"wrong" or unless another competing systematic effect con- 
spires to hide a bias this large. That the intrinsic scatter is in- 
deed larger than zero is, e.g., indicated by the observed seg- 
regation of cooling core clusters towards the high luminosity 
side for a given mass in theLx-M tot relation (Y. Chen et al., 
in prep.: lReiprich & Sar azin 2003) and by temporary luminos- 
ity increases in maj or mergers predicted by simulations (e.g., 
Rand all et al.l2002l) . We are currently studying these effects in 
more detail with better data. 

To derive an upper limit on the bias allowed by the WMAP 
3rd year data we used the largest £2 m and <x 8 values consis- 
tent with the q uoted uncertainties: £2 m = 0.27 and <x 8 = 0.79 
( Sperg eTet all200 6.note that we combined the uncertainties in 
Q. m h 2 and h in a worst case fashion, to be conservative). Then 
the quasi-mass function fit procedure was repeated for varying 
bias factors. The largest bias factor for which the 1-cr error el- 
lipse is still consistent with Q m = 0.27 and <x 8 = 0.79 is 1.5. 
Therefore, we conclude that the bias factor < 1.5, unless Q m 
and/or <x 8 are larger than allowed by the WMAP 3rd year data. 
Without being less conservative, a tighter constraint on the bias 
parameter seems currently not feasible with this method be- 
cause of the existing uncertainties in Q m and cr% and the strong 
sensitivity of the cluster mass function on these parameters. 

Despite the very good agreement with the new WMAP 
results, it should be noted that in their ROSAT and ASCA 
analysis RB02 had to use rather simplifying assumptions in 
their determination of cluster masses. For instance, the intr- 
acluster medium was assumed to be isothermal throughout 
the entire cluster volume and for some clusters only a rela- 
tively small number of photons was available for the determi- 
nation of the gas density profile. With the new generation of 
X-ray satellites these limitations can be overcome. Therefore, 
we are currently following-up all HIFLUGCS clusters with 
Chandra (Hudson et al., in prep.; Reipr icrTet al J 120061) and 
XMM-Newton (Nenestyan et al., in prep.). This should allow 
us to get a realistic estimate of the intrinsic scatter directly. 
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3. Summary of Conclusions 

1 . The shift in best fit values for £2 m and cr 8 between WMAP 
1st year and WMAP 3rd year data has a significant effect 
on predictions of cluster abundances because they are very 
sensitive to even small changes in these parameters. 

2. The WMAP 3rd year data are now consistent with the low 
cluster normalization found previously. 

3. The observed cluster X-ray luminosity-gravitational mass 
relation and X-ray luminosity function are in perfect agree- 
ment with mass functions predicted for Q m = 0.238 and 
cr 8 = 0.74, as favored by the WMAP 3rd year data. 

4. Assuming Q m = 0.238 and erg = 0.74, we find that the 
expected bias in normalization when determining the clus- 
ter luminosity-mass relation from flux-limited samples is 
small - much smaller than the factor of 2 claimed previ- 
ously. 

5. Assuming Q m = 0.27 and erg = 0.79, the largest values 
allowed by the new WMAP data, we derive a conservative 
upper limit on the bias factor of 1.5. 

6. A small bias indicates that the intrinsic scatter in the 
luminosity-mass relation is small. 

T he RB02 results, together with results from other w orks 
(e.g.. lBorgani et al.ll200lHseliakll2002t IViana et alJEool) that 
yielded the low cluster normalization, triggered a cr 8 -debate, 
which seems to have resulted in a loss of trust in cosmological 
constraints derived from clusters within the cosmology com- 
munity. Now that other well-respected methods like those that 
rely on measurements of anisotropics in the cosmic microwave 
background also tend to find lower values of <x 8 , we hope that 
trust in clusters will be restored. Clusters are ideally suited 
to address th e new fundamental questions in cosmology (e.g., 
lMohrll2005l) : planned X-ray cluster surveys like those to be 
performed with eROSITA 2 are expected to yield up to 100,000 
clusters out to large redshifts. The resulting detailed measure- 
ments of the evolution of cosmic structure are among the most 
promising routes to pin down the nature and evolution of dark 
energy. 
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